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Abstract 

We have studied the properties of the non-classical behavior of atoms in a 
double-slit interferometer. An indication of this behavior for metastable helium 
was reported by Kurtsiefer, Pfau and Mlynek [Nature 386, 150 (1997)] show- 
ing distinctive negative values of the Wigner function, which was reconstructed 
from the measured diffraction data. Our approach to explain this non-classical 
behavior is based on the de Broglie-Bohm-Vigier-Selleri understanding of the 
wave-particle duality and compatible statistical interpretation of the atomic wave 
function. It follows from the results that the atomic motion is non-classical be- 
cause it does not obey the laws of classical mechanics. However, there is no 
evidence that this atomic behavior violates the classical probability law of the 
addition of probabilities. 



PACS number: 03.65.Bz, 03.75.Dg, 03.75.* 

Key words: atomic interference, compatible statistical interpretation, (non)violation 
of the classical probability laws, Wigner's function 

Corresponding author: Mirjana Bozic, Institute of Physics, P.O.Box 57, 
11000 Belgrade, Yugoslavia, Tel: 381-11-3160260; Fax: 381-11-3162190; e-mail: 
bozicOphy . bg . ac . yu 



1 



I. Introduction 

The wave function ift(x, t) of the transverse motion of an atom in an atom inter- 
ferometer is a linear superposition of states with maxima at two spatially separated lo- 
cations. These states lead to negative values in Wigner's function W(x, p x , t) which is 
the quasi-probability distribution of coordinate x and momentum p x . The negative val- 
ues of W(x, p x , 0), reconstructed from measured and evaluated space distribution, were 
interpreted as a signature of the highly non-classical behavior of atoms in the atom inter- 
ferometer [1-3]. 

In this paper we study the properties and the cause of this non-classical behavior, 
using the compatible statistical interpretation (CSI) of a wave function [4-6]. The aim 
of this study is to clarify the meaning of the notion "non-classical motion (behavior)" of 
atoms in a double-slit interferometer. In our opinion, it is necessary to distinguish clearly 
two aspects of the notion "non-classical motion (behavior)". It may denote a motion 
(behavior) which does not obey the laws of classical mechanics and/or a motion which 
does not obey the classical probability laws, in particular the classical law of the addition 
of probabilities. 

We use CSI because the wave and corpuscular features of a wave function are incor- 
porated in a consistent manner into the basic statistical quantity of CSI, which is the de 
Broglian probability density P(x, p x , t). P(x, p x , t) is the probability density for a par- 
ticle, which is in the quantum state ift(x, t), to have a momentum p x and to be at x at 
the time t [5]. P(x, p x , t) satisfies both marginal conditions imposed by Wigner upon any 
joint probability distribution in phase space and it is always positive [6]. 

The coherence and the characteristic modulation of the momentum distribution found 
by Kaiser et al. [7] at the exit of a neutron interferometer was explained by Bozic and Marie 
[5], based on the P(x,p x ,t) function. Bozic and Arsenovic [10] compared the explanation 
of the same effect, based on Wigner's function and given by Lerner, Rauch, and Suda [8] 
and Suda [9], with an explanation based on the de Broglian probability density [5]. In 
this paper we use an analogous comparison of the time dependent Wigner's function to 
the de Broglian probability density. These are two different distributions in phase space, 
associated with the same wave function of the atomic transverse degree of freedom in the 
interferometer. 

In the following sections of this paper we summarize quantum properties of atomic 
motion in an interferometer. In Sec. II is written the solution of Schrodinger's equation 
for an interferometer in the form of the Fresnel-Kirchhoff diffraction integral, while in Sec. 
Ill is derived the time dependent wave function ip(x, t) of the transverse motion and, for 
a chosen set of parameters, graphs of the function \ip(x, i)\ 2 are presented. The transverse 
momentum distribution in the state ip(x, t) is evaluated and presented graphically in Sec. 
IV. The de Broglian probability density and Wigner's function in the state ip(x, t) are 
evaluated and graphically presented in Sects. V and VI, respectively. The comparison 
between these two probability distributions is also given in Sec. VI, while in Sec. VII are 
derived concluding remarks about the properties of non-classical behavior of atoms in an 
interferometer. 
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II. The application of the Fresnel-Kirchhoff diffraction formula 



We want to determine the wave function of the transverse motion of an atom which 
travels with velocity v = vi y through region I (see Fig. 1), towards the slits and is then 
sent through the slits to region II. For this reason we shall determine in this section a 
stationary solution of the time-dependent two dimensional Schrodinger equation 

h 2 ( d 2 d 2 \ d 

+ —)y(x,y,t)=in—y(x,y,t). (1) 



2m \dx 2 dy 2 J ' ' dt 
The stationary solution of Eq. (1) has the form 

y(x,y,t) = e- i " t $(x,y), (2) 

where fiu = p 2 /2m and p = mv = Tik. The space dependent function 3>(x, y) satisfies the 
equation 

+ — T J y) = hu$(x, y). (3) 



2m \dx 2 dy 2 
The solution of Eq. (3) in region I is a spherical wave 

Apikr' 

$(P') = $( X >,y>) = —^, (4) 

where A is a constant and r' is the distance (Fig. 1) from the source (Pq) to the point P' = 
(x', y') in region I. The spherical wave at the slit points (x', y' = 0) may be approximated 
by a plane wave, since the distance a of the double-slit screen from the source Po is very 
large compared to the width of the slits. Consequently, without a loss of generality, for 
$(x', y' = 0) at the border of region I we may choose the function 

for one open slit, and the function 

-%>x'>-±-6 
Mx',0) = { 1/y /2S t a +6 > x ,>a (56) 

all other values of x' 

for two open slits. This means that in region II the solution of Eq. (3) is given by the 
formula of the Fresnel-Kirchhoff diffraction [11] 

\A e lka f e iks 

Hx,y) = -7rc [1 + cosx], (6) 

za a j s 

where s = ^/y 2 + (x f — x) 2 , cosx = y/s, A = 2n/k, while A = {x'; —(A/2) — 5 < x' < 
— (A/2)} when the lower slit is open and upper slit is closed, and A = {x'\ (A/2) < x' < 
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(A/2) + 5 or -(A/2) - 5 < x' < -(A/2)} when the two slits are open. The constant A 
will be chosen from the normalization condition. 

The spatial distribution of the transverse degree of freedom as a function of evolution 
time was investigated in a double slit experiment [1,2] with metastable helium atoms. We 
will apply the Fresnel-Kirchhoff diffraction formula to analyze the experiment of Kurtsiefer, 
Pfau, and Mlynek [1]. A diagram of the apparatus is shown in Fig. 2. Atoms are emitted 
from a gas-discharge source operating in the pulse operation mode. The beam is collimated 
by a 5 — /um-wide slit and then is sent through a double-slit structure with a slit separation 
A + 5 = 8 fim and an opening 6 = 1 fxm. The atoms then propagate for a distance d to a 
time- and space-resolving detector. Atom beam velocities lie between 1000 and 3000 ms" 1 . 
We shall use the parameters of this experimental arrangement for the following calculations. 

III. Time-dependent wave function of the transverse motion 

Assuming that the motion of an atom along the y-axis can be treated classically and 
that the transverse motion is quantized, one may use the relation y = vt and determine 
the time dependent function of the transverse motion from the function $(x, y), by the 
following definition 

$(x,y) = Q(x,vt) = i/>(x,t). (7) 

The graphs of the function |<&(x, y)\ 2 = \ip(x, t)\ 2 for k = An ■ 10 10 m -1 and for the chosen 
set of values of the coordinate y (t = my/Tik) are presented in Figs. 3 and 4. 

Very close to the slit on the single-slit graphs (Fig. 3) we see the minima of the wave 
function for x = x c , where x c = —4 fim is the coordinate of the slit center. But, with 
increasing y, the maximum is present at x = x c for all y. This maximum becomes wider 
and wider with increasing y. 

On double-slit graphs (Fig. 4) we clearly see that near the slits the wave function 
consists of two widely separated Gaussians on which small oscillations are superimposed. 
With increasing distance from the slits the Gaussian-like maxima spread and start to 
overlap, so that the third maximum with superimposed oscillations start to develop. This 
region of y corresponds to Fresnel diffraction. With further increase of y (i), distinct 
equally spaced oscillations develop, which correspond to the Fraunhofer diffraction limit. 

IV. The transverse-momentum distribution 

The time dependent function defined by Eq. (7) should be a solution of the one-di- 
mensional time-dependent Schrodinger's equation. Therefore, we may assume that it can 
be written in the form 

-I />oo -I poo 

iKM) = ^=== / c(p x )e^ h e-^ t dp x = -= / c'ik^e^e-^dk^ (8) 

V2"7m J-oo V27T J-ao 

where \c(p x )\ 2 dp x = \c'(k x )\ 2 dk x = 1, p x = hk x , c'(k x ) = Vhc(p x ) and Tiu x = 
p 2 /2m. From Eq. (8) we may determine the transverse-momentum distribution |c(p x )| 2 = 
\c'(k x )\ 2 /h in the state ip(x,t). At first, one determines 

1 f°° 

C{k x ,t) = ^= i;(x,t)e-^ x dx (9) 
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by performing the Fourier-transform of the function ip(x,t), defined by Eq. (7), taking t 
as a parameter. If Eq. (8) is valid, then it should be 

C(k x ,t)=c'(k x )e-^ t . (10) 

Consequently, 

\c'(k x )\ 2 = \C(k x ,t)\ 2 . (11) 

The graph of \c'(k x )\ 2 = fi\c(p x )\ 2 for one slit is given in Fig. 5a and for two slits in Fig. 
5b. 

Our numerical calculation for various values of t, show that \c'(k x )\ 2 is independent 
of t. This fact justifies the assumptions of Eqs. (7) and (8) as well as the statement of 
Kurtsiefer, Pfau, and Mlynek [1] that the longitudinal motion of the atoms at velocities v 
of several thousand meters per second can be treated completely classically. 

We compared also the transverse momentum distribution ^'(/c^)! 2 (evaluated as de- 
scribed above and presented in Fig. 4) with the absolute value square of the Fourier trans- 
form 

Fi(k x ) = -= \ ( p i (x\0)e-^ x 'dx' (12) 

V27T J-oo 

of the function <j>i(x', 0), i = 1, 2. After the evaluation of the latter integral one finds 



le 



ik x A/2 



and 



„ n s 2 . k x S k x (A + d) 

F 2{k x ) = - — 1= sin — cos : 

k x yno * * 



I t-i / t \ 1 2 4 -2 k xS 2 k x (A + S) 

|F2(MI = khTs sm — cos —2~- (14) 



We found that \c'(k x )\ 2 for one slit is practically identical to \Fi(k x )\ 2 and that \c'(k x )\ 2 
for two slits is practically identical to \F2(k x ) \ 2 ■ 

By comparing the spatial distributions for one and two slits shown in Figs. 3 and 
4, one must conclude that the presence of the second slit influences the motion of each 
atom, independent of the slit through which it has passed to region II (see Fig. 1). This 
influence is also very well seen by comparing the momentum distributions for one and two 
slits, presented in Fig. 5. Certain values of the particle's transverse momentum, which 
are allowed with one slit, are not allowed when both slits are open. This fact is also a 
signature of a non-classical atomic behavior that can be understood in a similar way to the 
quantization of the electronic orbits in atom based on de Broglie's wavelength. It appears 
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that the atomic matter wave excludes certain values of transverse momentum and favors 
others, which is an evident quantum effect. 

V. The De Broglian probability density 

According to the CSI of a wave function, in an ensemble of particles in a pure state 
presented by Eq. (8), different particles may have different momenta. Recall that the 
probability density of p x is (c^)! 2 . However, each particle is surrounded by the same 
wave [5]. In other words, a particle and a wave are two different, but compatible, entities. 

The de Broglian probability density, P(x,p x ,t), of a particle in the quantum state 
ip(x, t) is the probability density for the particle to have a momentum p x and to be at a 
position x at time t [5,10], 

P(x,p x ,t) = \^{x,t)\ 2 \c{p x )\ 2 = P'(x,k x ,t)/h = \i;(x,t)\ 2 \c'(k x )\ 2 /h. (15) 

P(x,p x ,t) satisfies both marginal conditions 

t)dp x = \ip(x,t)\ 2 , (16) 

j P(x,p x ,t)dx = \c(p x )\ 2 (17) 

imposed by Wigner upon a joint probability distribution in phase space [12]. For operators 
having the form F(x,p x ) = F\{x) + Fz(p x ), the probability density P(x,p x ,t) satisfies 
Wigner's condition that the quantum mechanical average value of an operator is equal to 
the classical average value of the corresponding classical function. 

Both P(x, p x , t) and Wigner's function are determined by the state ip(x, t). Unlike the 
Wigner function, P(x,p x ,t) is always positive. Despite this fact, P(x,p x ,t) also reflects 
the non-classical behavior of atoms in the state ip(x,t). 

Since the simultaneous measurement of a coordinate and momentum is not possible, 
P(x,p x ,t) cannot be measured in a single experiment. However, one could experimentally 
determine the probability density of a coordinate x and momentum p x in the state ip(x, t), 
i.e. P(x,p x ,t), by measuring separately the distributions \ift(x,t)\ and |c(p x )| 2 . These 
distributions reflect the non-classical behavior, as pointed out in the previous section. 

VI. The De Broglian probability distribution and Wigner's function 

In Figs. 6 and 7 we present the graphs of the de Broglian probability density of 
a coordinate x and transverse momentum p x , P(x, p x , t), for y = 120 mm (t = y/v = 
6.01 x 10" 5 s) and y = 240 mm (t = y/v = 12.02 x 10 5 s). For the same values of y we 
present in Figs. 8 and 9 the plots of the Wigner distribution function, evaluated from the 
definition [12,13] expression 

W(x, p x , *) = -?-/ dxe 2ip ^ /n ij*(x + x, t)ip(x -x,t) = W'(x, k x , t)/h. (18) 
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It is clear from Figs. 6-9 that W(x, p x , t) and P(x, p x , t) are very different functions. 
Consequently, from their forms and properties are derived different interpretations of the 
behavior of quantum particles. It was shown by Janicke and Wilkens [14] and Kurtsiefer, 
Pfau, and Mlynek [1] that Wigner's function W(x, p x , 0) may be reconstructed from eval- 
uated and measured values of \ip(x, t)\ 2 for various values of t. The negative values of 
W(x, p x , 0) were interpreted as a signature of an atom's non-classical behavior. These 
negative values are also associated with the requirement of Heisenberg's uncertainty rela- 
tionship that a quantum mechanical particle has to be described by an area of uncertainty 
in phase space no smaller than AxAp x = K/2 [3]. The authors pointed out that the 
negative values reflect the impossibility of joint measurement of position and momentum. 

However, we interpret de Broglian probability density, presented in Figs. 6 and 7, as 
an objective probability density of particle coordinate and momentum. The impossibility 
of simultaneous measurements of a particle's x and p x does not forbid us from assuming 
that their joint distribution objectively exists. The important fact is that this assumption 
does not lead to any contradiction with the facts derived from measurable distributions. 
One can see that this joint probability density is consistent with the measurable probability 
density of position and the measurable probability density of momentum. For example, 
for values of p x for which |c(p x )| 2 = 0, the joint distribution P(x, p x , t) is also equal to 
zero. Thus, if there is no particle with a certain value of momentum p x , this value can not 
be found anywhere during the measurement of momentum. Similar reasoning is valid for 
space points x in which \ip(x, t)\ 2 = 0, since P(x,p x , t) is also equal to zero in these space 
points for any value of momentum. Therefore, at a point x no particle will be detected in 
the experiment. 

One can see in Figs. 8 and 9 that Wigner's function W(x,p x ,t) may take values 
different from zero at the points x and p in which either ip(x, t) = or c(p x ) = 0. Despite 
this property, inconsistent with a notion of a joint probability, the Wigner function satisfies 
the marginal conditions stated by Eqs. (16) and (17). It is well known that Wigner's 
function may assume negative values, even though it is a joint probability distribution by 
definition. Because of this, it is possible to satisfy Eqs. (16) and (17). Thus, two different 
properties of Wigner's function, inconsistent with a notion of a joint probability, cancel 
each other and make it possible to satisfy two marginal conditions. This is clearly seen by 
comparing results presented in Figs. 6,7 and 8,9. We note in Fig. 8,9 the negative peaks 
in the x-dependence of the Wigner function for those values p x of momentum for which 
|cfe)| 2 = 0. 

VII. Conclusion 

The properties and cause of non-classical behavior of atoms in the atomic interferom- 
eter are studied using the stationary solution $(x, y) of the two-dimensional Schrodinger's 
equation. The solution was written in the form of the Fresnel-Kirchhoff diffraction integral. 
The time dependent wave function of the transverse motion was derived from the Fresnel- 
Kirchhoff diffraction integral, using the relation $>(x, y) = $>(x, vt) = ip(x, t), where v is 
the initial longitudinal atomic velocity. The latter relation was used in Refs. [1,2], where 
it was justified by the experimental facts, suggesting that the longitudinal atomic motion 
was classical and that the transverse motion was quantum. We determined the transverse 
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momentum distribution in the state ip(x, t), by evaluating its Fourier transform. 

We calculated \ip(x, t)\ 2 for one and two-slit interferometers, and presented results in 
Figs. 3 and 4. From the data one can see that the evolution (and spreading) of waves from 
different single slits and their interference (overlap) determine Fresnel's and Fraunhofer's 
regimes and the transition from the former to the latter. We conclude that this spatial 
distribution, which reflects the non-classical atomic motion, is due to a real atomic wave 
that is associated with each atom and that influences its motion. The obstacle in front of 
the incoming atoms determines the concrete form of this influence. 

Our results show that the de Broglie [15], the Bohm and Vigier [16] and the Selleri 
[17] understandings of wave-particle duality is applicable to the explanation of the non- 
classical motion of atoms in an interferometer. The conclusion that the motion is non- 
classical means that it is different from the motion of a particle which obeys the laws 
of classical mechanics. This difference is due to the fact that with a classical particle 
no wave is associated, whereas the atom is accompanied by its wave. This conclusion 
follows from the measured atomic distribution [1,2], and its theoretical explanation, based 
on the particle's wave function, in this and in the previous works [1,2,18]. Therefore, 
the application of methods for determination of the amplitude and phase structure of the 
atomic wave field, like the method of Raymer, Beck, and McAlister [19] would be of great 
importance. 

However, neither from the measured space distribution [1,2], nor from the transverse 
momentum distribution evaluated in this paper, does it follow that the motion of atoms 
in the double-slit interferometer violates the classical probability laws, especially the law 
of the addition of probabilities. From the fact that Wigner's function, associated with the 
state ip(x, t), takes a negative value if does not follow that atomic motion violates the latter 
law. It rather seems that the classical probability laws are satisfied and are consistent with 
the de Broglie-Bohm-Vigier-Selleri [15-17] understanding of the wave-particle duality and 
the compatible statistical interpretation of a wave function [4-6]. 

To the best of our knowledge, the atomic transverse momentum distribution in the 
interferometer has not been measured. However, as pointed out in this paper, it is very 
important characteristic of the quantum state. Experimental evidence of the transverse 
momentum distribution would contribute a lot to our understanding of the quantum nature 
of atomic motion. 
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Figure captions 



Figure 1. Illustration of the diffraction formula presented with Eq. (6). 

Figure 2. Diagram of apparatus used in Ref. [1] to observe atomic interference patterns. 

Figure 3. The function \ip(x, t)\ 2 = |<&(x, y = vt)\ 2 for a single slit evaluated from 
Eq. (6) with A chosen, from the condition J \<&(x, y = vt)\ 2 dx = 1 for a given y. Other 
parameters are: k = 4vr • lO^m -1 , v = hk/m = 1995.58 m/s, m = 6.64632 • 10" 27 kg is 
the mass of the Helium atom. 

Figure 4. The function \ip(x, t)\ 2 = \$>(x, y = vt)\ 2 for a double-slit evaluated from 
Eq. (6) with A chosen, from the condition J \<&(x, y = vt)\ 2 dx = 1, for a given y. Other 
parameters are: k = 4tt ■ 10 10 m -1 , v = hk/m = 1995.58 m/s, m = 6.64632 ■ 10 -27 kg is 
the mass of the Helium atom. 

Figure 5. Momentum distribution \c'(k x )\ 2 = \c(p x )\ 2 ■ h in the state ip(x, t) with 
parameters given in the caption of Figs. 3 and 4. a) One open slit, b) Two open slits. 

Figure 6. The de Broglian probability density P'(x, k x , t) = hP(x, p x , t) in the single 
slit state ip(x, t) with parameters given in the caption of Fig. 3. 

Figure 7. The de Broglian probability density P'{x, k x , t) = %P{x, p x , t) in the 
double-slit state ip(x, t) with parameters given in the caption of Fig. 4. 

Figure 8. Wigner's function W'(x, k x , t) associated with the single slit state ip(x, t) 
and evaluated from (18). Parameters are given in the caption of Fig. 3. 

Figure 9. Wigner's function W'(x, k x , t) associated with the double-slit state ip(x, t) 
and evaluated from (18). Parameters are given in the caption of Fig. 4. 
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